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This paper reviews the constraints imposed on the solar neutrino mixing parameters by data collected by the 
Sudbury Neutrino Observatory (SNO). The SNO multivariate analysis is reviewed. The global solar neutrino 
analysis is emphasized in terms of matter-enhanced oscillation of two active flavors. An outline of how SNO 
uses the data to produce oscillation contour plots and how to include the relevant correlations for the new salt 
data in similar oscillation analyses is summarized. 
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1. Introduction 

The deficit of detected neutrinos coming from the 
Sun compared with our expectations based on lab- 
oratory measurements, known as the Solar Neutrino 
Problem, was one of the outstanding problems in basic 
physics for over thirty years. It appeared inescapable 
that either our understanding of the energy producing 
processes in the Sun was seriously defective, or neu- 
trinos, one of the fundamental particles in the Stan- 
dard Model, had important properties which had not 
been measured. It was indeed argued by some that 
we needed to change our ideas on how energy was 
produced in fusion reactions inside the Sun. Others 
suggested that the problem arose due to peculiar char- 
acteristics of neutrinos such as vacuum or matter os- 
cillations. It is useful to review the evolution of our 
understanding from the data collected by various so- 
lar neutrino experiments. The new analysis of the 
salt data collected by the Sudbury Neutrino Obser- 
vatory (SNO) [l| will be described, together with the 
technique used to combine the results of many solar 
neutrino experiments. 



2. Solar Neutrinos 

The energy in the Sun is produced by nuclear reac- 
tions that transform hydrogen into helium. Through 
the fusion reactions, four protons combine to form a 
helium nucleus containing two protons and two neu- 
trons. The only reactions that allow this to happen 
are caused by weak interactions like nuclear beta de- 
cay. Each time a neutron is formed, there must be an 
associated positron and electron neutrino produced. 
Neutrinos can travel directly from the core of the Sun 
to the Earth in a about eight minutes and hence pro- 
vide a direct way to study thermonuclear processes in 
the Sun. The detailed predictions of the solar electron 
neutrino flux have been produced by John Bahcall and 
his collaborators from the 1960's until now. Their cal- 
culations are refereed to as the Standard Solar Model 
(SSM). In this proceeding, the Bahcall-Pinsonneault 
calculations are compared to experimental results. 



It is known that neutrinos exist in different fla- 
vors corresponding to the three charged leptons: the 
electron, muon, and tau particles. If neutrinos have 
masses, flavor can mix and a neutrino emitted in 
a weak interaction is represented as a superposition 
of mass eigenstates. In the case of three flavors of 
neutrino, the mixing matrix U is called the Maki- 
Nakagawa-Sakata-Pontecorvo (MNSP) matrix and 
Vi = YliUti\vi)- Here the neutrino mass eigenstates 
are denoted by V\ with % = 1,2,3, while the flavor 
eigenstates are labeled (e,/i, r). The most general 
form of mixing for three families of neutrinos can be 
simplified so that only two neutrinos participate in the 
oscillations. Hence, the survival probability for solar 
neutrinos propagating in time takes the approximate 
form 
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The mixing angle is represented by 6, L is the distance 
between the production point of u e and the point of 
detection of up, E is the energy of the neutrino, and 
Am 2 = m 2 — m 2 is the difference in the squares of the 
masses of the two states Uj and v\ which are mixing. 
The function 5 e p is the usual Kronecker delta. The 
numerical constant 1.27 is valid for L in meters, E in 
MeV, and Am 2 in eV 2 . The energy of a neutrino de- 
pends on the type of nuclear reaction which produced 
it. By studying the evolution of the solar neutrinos as 
a function of L, all the physics is embedded in one an- 
gle 9, one mass difference Am 2 , and the sign of Am 2 . 
This corresponds to the extraction of the three MNSP 
elements: U e \, U e 2, and U e ^. 



3. Sudbury Neutrino Observatory 

The Sudbury Neutrino Observatory (SNO) is a 
1,000 ton heavy-water Cerenkov detectorjj] situated 
2 km underground in INCO's Creighton mine in 
Canada. Another 7,000 tons of ultra-pure light water 
is used for support and shielding. The heavy water is 
in an acrylic vessel (12 m diameter and 5 cm thick) 
viewed by 9,456 PMT mounted on a geodesic structure 
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18 m in diameter; all contained within a polyurethane- 
coated barrel-shaped cavity (22 m diameter by 34 m 
high). The solar-neutrino detectors in operation prior 
to SNO were mainly sensitive to the electron neutrino 
type; while the use of heavy water by SNO allows 
neutrinos to interact through charged-current (CC), 
elastic-scattering (ES), or neutral-current (NC) inter- 
actions. The determination of these reaction rates is 
a critical measurement in determining if neutrinos os- 
cillate in transit between the core of the Sun and their 
observation on Earth. 

During the pure D2O phase of the experiment, the 
signal was determined with a statistical analysis based 
on the direction, cos# sun , the position, R, and the ki- 
netic energy, T, of the reconstructed events assuming 
the SSM energy spectrum shape The final se- 
lection criteria were T > 5 MeV and R < 550 cm. 
The result of the extended maximum-likelihood fit 
yields ja] 

$cc = 1.76lgS18S xlO'cm-V 1 , 

$ ES - 2.39 ±8:it±8:ii xio'cm-V 1 , (2) 

<I>nc = 5.09 xlO'cm-V 1 . 

The excess of the NC flux over the CC and ES fluxes 
implies neutrino flavor transformations. There is also 
a good agreement between the SNO NC flux and the 
total 8 B flux of 5.05±o:gi xlO 6 cm^s" 1 predicted by 
the SSM. A simple change of variables that resolves 
the data directly into electron and non-electron com- 
ponents [6j indicates clear evidence of solar neutrino 
flavor transformation at 5.3 standard deviations 

e = 1.76t° : °ft°;» Xl0 6 cm-V\ (3) 
V - 3.41±g;«±g : « xlO^m-V 1 . (4) 

Allowing a time variation of the total flux of solar 
neutrinos leads to day/night measurements by SNO, 
which are sensitive to the neutrino type Q 

A DN (total) = (-24.2 ± 16.1 ±H) % , (5) 
A DN (e) = (12.8 ± 6.2 t\i)%. (6) 

By forcing no asymmetry in the <f> e + 4>^ T rate, i.e. 
Aqn (total) = 0, the day/night asymmetry for the 
electron neutrino is A DN (e) = (7.0 ± 4.9 tx.f)- 

SNO published its first results of the salt phase 
in coincidence with the PHYSTAT2003 conference. 
The measurements were made with dissolved NaCl in 
the heavy water to enhance the sensitivity and signa- 
ture for neutral-current interactions. Neutron capture 
on CI typically produces multiple 7 rays while the 
CC and ES reactions produce single electrons. The 
greater isotropy of the Cerenkov light from neutron 
capture events relative to CC and ES events allows 
good statistical separation of the event types. The 
degree of the Cerenkov light isotropy is determined 
by the pattern of PMT hits. This separation allows 



a precise measurement of the NC flux to be made in- 
dependent of assumptions about the CC and ES en- 
ergy spectra. To minimize the possibility of introduc- 
ing biases, SNO performed a blind analysis for the 
model independent determination of the total active 
8 B solar neutrino. In this analysis, events are statisti- 
cally separated into CC, NC, ES, and external-source 
neutrons using an extended maximum- likelihood tech- 
nique based on the distributions of isotropy, cos9 sun , 
and radius, R, within the detector. To take into 
account correlations between isotropy and energy, a 
2D joint probability density function (PDF) is con- 
structed. This analysis differs from the analyses of the 
pure D2O data [(JlI since (1) correlations are explic- 
itly incorporated in the signal extraction and (2) the 
spectral distributions of the ES and CC events are not 
constrained to the S B shape, but are extracted from 
the data. Cerenkov event backgrounds from /3 — 7 
decays are reduced with an effective electron kinetic 
energy threshold T > 5.5 MeV and a fiducial volume 
with radius R < 550 cm. 

The extended maximum-likelihood analysis gives 
the following & B fluxes [lj 

<T> — 1 +0.08 +0.06 v in6 

*CC — 1-59-0.07 -0.08 x 10 cm s , 

$es = 2.21 toil ± 0.10 x 10 6 cm-V 1 , (7) 
$nc = 5.21 ± 0.27 ± 0.38 x 10 6 cm" 2 s _1 . 

The systematic uncertainties on the derived fluxes are 
shown in Tabled These fluxes are in agreement with 
previous SNO measurements and the SSM. The ra- 
tio of the 8 B flux measured with the CC and NC re- 
actions then provides confirmation of solar neutrino 
oscillations 

= 0.306 ± 0.026 ± 0.024 . (8) 

S'NC 



4. How to Use the SNO Data 

The SNO CC, ES and NC fluxes are statistically 
correlated, since they are derived from a fit to a sin- 
gle data set. The statistical correlation coefficients 
between the fluxes in the salt phase are 

Pccnc = -0.521, 

Pcc.es = -0.156, (9) 

Pes.nc = -0.064. 

These can be used with the statistical uncertainties 
quoted by SNO Q to write down the statistical co- 
variance matrix for the salt fluxes. Systematic un- 
certainties between fluxes can be correlated as well. 
Some sources of systematic error, such as neutron 
capture efficiency, affect only one of the three fluxes, 
and so can be considered to be uncorrelated with the 
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Source 


NC 


CC 


ES 


Energy scale 


-3. 7, +3. 6 


-1.0, ±1.1 


±1.8 


Energy resolution 


±1.2 


±0.1 


±0.3 


Energy non-linearity 


±0.0 


-0.0, ±0.1 


±0.0 


Radial accuracy 


-3.0, ±3. 5 


-2. 6, ±2. 5 


-2. 6, ±2. 9 


Vertex resolution 


±0.2 


±0.0 


±0.2 


Angular resolution 


±0.2 


±0.2 


±2.4 


Isotropy mean 


-3. 4, ±3.1 


-3. 4, ±2. 6 


-0.9, +1.1 


Isotropy resolution 


±0.6 


±0.4 


±0.2 


Radial energy bias 


-2. 4, ±1.9 


±0.7 


-1.3, ±1.2 


Vertex Z accuracy 


-0.2, ±0.3 


±0.1 


±0.1 


Internal neutrons 


-1.9, ±1.8 


±0.0 


±0.0 


Internal background 


±0.1 


±0.1 


±0.0 


Neutron capture 


-2.5,±2.7 


±0.0 


±0.0 


Cerenkov backgrounds 


-1.1,±0.0 


-1.1,+Q.O 


±0.0 


AV events 


-0.4,±0.0 


-0.4,±0.0 


±0.0 


Total uncertainty 


-7.3,±7.2 


-4.6,±3.8 


-4.3,±4.5 



Table I Systematic uncertainties (in %) on fluxes for the 
spectral shape unconstrained analysis of the salt data set. 



other fluxes. Other systematics can be either 100% 
correlated (e.g. radial accuracy) or 100% anticorre- 
lated (e.g. isotropy mean). The most important anti- 
correlated systematic is the isotropy mean. Isotropy 
is important for separating CC and ES events from 
NC events, so CC and ES will have a negative cor- 
relation with the NC flux (and a positive correlation 
with each other) for the isotropy uncertainty. Tabic ITT1 
shows the sign of the correlation for each systematic of 
Table [I] Using the table of systematics and the signs 
for the correlations, one can assemble an individual 
covariance matrix for each systematic. Then, to get 
the total covariance matrix for the CC, ES and NC 
fluxes, one simply adds all of the covariance matrices 
together. 

Even when fluxes are being analyzed as opposed to 
energy spectra, it is best to determine the effect of 
energy-related systematics at each grid point in the 
Am 2 — tan 2 plane. For the salt analysis, these in- 
clude energy scale and energy resolution; the uncer- 
tainty due to energy non-linearity is tiny so that it can 
reasonably be ignored. The energy scale uncertainty 
is implemented as a 1.1% uncertainty in the total en- 
ergy; while the energy resolution has an uncertainty 
which is energy dependent for T > 4.975 MeV 



Aery 

0~t 



0.035 ± 0.00471 x (T — 4.975) 



(10) 



and = 0.034 for T < 4.975 MeV. Here T is the 
reconstructed kinetic energy. For all other systemat- 
ics, it is assumed that the effect on the fluxes is the 
same for all oscillation parameters. 

When SNO quotes $cc = 1-59 x 10 6 cm^s- 1 , it 
refers to the integral flux from zero to the endpoint as- 



Source 


NC 


CC 


ES 


Energy scale 


-1-1 


-1-1 


+1 


Energy resolution 


+ 1 


+ 1 


+ 1 


Energy non-linearity 


+1 


+1 


+ 1 


Radial accuracy 


+1 


+1 


+1 


Vertex resolution 


+ 1 


+ 1 


+1 


Angular resolution 


+1 


+1 


-1 


Isotropy mean 


+1 


-1 


-1 


Isotropy resolution 


±1 


+1 


+ 1 


Radial energy bias 


+ 1 


+ 1 


+ 1 


Vertex X accuracy 


+1 


-1-1 


+1 


Vertex Y accuracy 


+1 


+1 


+1 


Vertex Z accuracy 


+1 


-1 


- 1 


Internal neutrons 


+ 1 


o 


o 


Internal background 


±1 


±1 


±1 


Neutron capture 


±1 








Cerenkov backgrounds 


±1 


±1 


±1 


AV events 


±1 


±1 


±1 



Table II Signs of systematic correlations, relative to its 
effect on the NC flux. An entry of ±1 indicates a 100% 
positive correlation, —la 100% negative correlation, and 
means no correlation. 



suming an undistorted 8 B spectrum. It implies that 
the number of events attributed to CC interactions 
above T — 5.5 MeV is equal to the number of events 
that would be observed if the v e flux follows the 8 B 
spectral shape. The S B spectral shape aspect of this 
definition is only for normalization; there is no as- 
sumption of any spectral shape when extracting the 
number of events during the salt phase. Similar defi- 
nitions apply for the NC and ES fluxes. 

For the comparison of the SNO CC rate with the 
theoretical rates for a set of oscillation parameters, 
the <I>cc nux is 



Sb \ <hsM(E v )dE v S(T,T e ,E v ), 
Jo 

with the scale S(T, T e ,E v ) is equal to 

Jo°° Jo" JTs F(T, T e , E y )P ee {E v )dTdT e dE y 
Jo°° Jo°° JTs F ( T > T ei E v )dTdT e dE„ 

where 



(11) 



(12) 



F(T, T e , E v ) = <fi SSM (E v ) da ^T e ) N{T ^ 4) 

(13) 

The factor fs allows the total 8 B solar neutrino flux 
to float from the SSM value, E v is the neutrino en- 
ergy, P ee is the survival probability, T e is the true 
recoil electron kinetic energy, and T is the observed 
electron kinetic energy; while N(T e ,a^) is a Gaus- 
sian energy response function for T with o~t(T) = 
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-0.145 + 0.392\/r + 0.0353T. It is a similar definition 
for the SNO ES flux, remembering to include the con- 
tribution from !A, T using the appropriate cross section 
and (1 — P ee )- There is no ambiguity in interpreting 
NC flux since it is equal to the total SSM flux. 



5. Global Fits 

This section summarizes the constraints from so- 
lar neutrino data in a global analysis. The allowed 
region in the oscillation Am 2 — tan 2 plane is ob- 
tained by comparing the measured rates to the calcu- 
lated SSM solar neutrino rate. We consider a set of 
N observables R n for n = 1,2, ■■■ ) N with the as- 
sociated set of experimental observations i?° xp and 
theoretical predictions i?^ h . In general, one wants to 
build a x 2 function which measures the differences 
i?° xp — i?^ in units of the total experimental and the- 
oretical uncertainties. This task is completely deter- 
mined from the estimated uncorrelated errors u n and 
a set of correlated systematic errors caused by K 
independent sources. The correlation coefficients be- 
tween the different observables are p(u n , u rn ) = ±5 nm 
and p(c^,c^) = ±5kh- The covariance matrix takes 

the form a 2 m — S nm u n u m + J2k=i c n c m aim au the 
experimental information is combined together in a 
global x 2 

N 

xL= E (c-OkJ^^r-O. (14) 

n,m — 1 

The salt shape-unconstrained fluxes presented 
here, combined with shape-constrained fluxes and 
day/night energy spectra from the pure D2O phase 0, 
0, place impressive constraints on the allowed neu- 
trino flavor mixing parameters. In the fit, the ratio 
fs of the total 8 B flux to the SSM value is a free pa- 
rameter together with the mixing parameters. A com- 
bined x 2 n t to SNO D2O and salt data alone yields 
the allowed regions in Am 2 and tan 2 9 shown in Fig.Q] 
There are certainly correlations between the salt and 
the D2O phase, since it's the same detector. However, 
these correlations are estimated to be negligibly small. 

The Xcov calculated above from the SNO NC, CC 
and ES fluxes is added to a global analysis which in- 
cludes data from all the other solar neutrino experi- 
ments. Systematic errors that are correlated between 
different experiments, such as cross section uncertain- 
ties or uncertainties on the 8 J5, are accounted for by 
including the covariance terms between different ex- 
perimental results. The effect of the 8 B spectral shape 
uncertainty is determined at each grid point in the os- 
cillation plane. 

The global analysis includes the Homestake re- 
sults 0, the updated Gallium flux measurements 
Hcj . the SK zenith spectra |ll|, and the D2O and salt 



,oio' 3 1 — — — — g 

> E E 




E 90% CL 

!0 f 95% CL 

jd : 99% CL 

10 T 99.73% CL 



]f , 1 , — 1 — 1 — 1 — 

iu -3 -2 -1 

10 10 10 1 10 

tan 2 6 

Figure 1: SNO-only neutrino oscillation contours, 
including pure D2O day/night spectra, salt CC, NC, ES 
fluxes, with 8 B flux free and hep flux fixed. The best-fit 
point is Am 2 = 4.7 x 1CT 5 , tan 2 9 = 0.43, f B = 1.03, 
with x 2 /d.o.f.=26.2/34. The inside of the covariance 
regions is allowed. 

results from SNO 0,0,0- At each grid point in the 
Am 2 — tan 2 9 plane, the expected rate for each energy 
bin is calculated and compared to the measured rate. 
The free parameters in the global fit are the total S B 
flux, the difference of the squared masses Am 2 , and 
the mixing angle 9. The higher energy hep v e flux 
is fixed at 9.3 x 10 3 cm' 2 s _1 . Contours are gener- 
ated in Am 2 and tan 2 9 for Ax 2 ov = 4.61 (90% CL), 
5.99 (95% CL), 9.21 (99% CL), and 11.83 (99.73% 
CL). We assume a Gaussian distribution of R^ p for 
a given value of the true parameters dm 2 and tan 2 9 
when we map the survival probability into the MSW 
plane [12j. As presented in FigEJa), the combined 
results of all solar neutrino experiments can be used 
to determine a unique region of the oscillation pa- 
rameters; the allowed region in this parameter space 
shrinks considerably to a portion of the Large Mixing 
Angle (LMA) region. 

A global analysis including the KamLAND reactor 
anti-neutrino results |13| shrinks the allowed region 
further, with a best-fit point of Am 2 = 7.lj~ Q ' 6 x 
10 -5 eV 2 and 9 = 32.5^2 3 degrees, where the er- 
rors reflect la constraints on the 2-dimensional re- 
gion. This is summarized in Fig.|2jb). With the new 
SNO measurements, the allowed region is constrained 
to only the lower band of LMA at > 99% CL. The 
best-fit point with a one dimensional projection of the 
uncertainties in the individual parameters (marginal- 
ized uncertainties) is Am 2 = 7.1^' 3 x 10~ 5 eV 2 and 
9 = 32.51^ de grees. This disfavors maximal mix- 
ing at a confidence level equivalent to 5.4 standard 
deviations and indicates tan 2 9 < 1 . In our interpre- 
tation, the xlov f° r ^ = 45.0 is 5.4 2 higher than the 
best LMA fit. The solution tan 2 9 < 1 corresponds to 
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tan 2 9 tan 2 9 



Figure 2: Allowed region of the Am 2 — tan 2 8 plane 
determined by a \ 2 fit to (a) the Chlorine, Gallium, SK, 
and SNO experiments. The best-fit point is 
Am 2 = 6.5 x 10~ 5 , tan 2 6 = 0.40, f B = 1.04, with 
X 2 /d.o.f.=70.2/81. (b) Solar global + KamLAND. The 
best-fit point is Am 2 = 7.1 x 10" 5 , tan 2 9 = 0.41, 
/b = 1.02. The inside of the covariance contours is the 
allowed region. 



the neutrino mass hierarchy m,2 > mi . 



ergy dependence of solar oscillations with Large Mix- 
ing Angle (LMA) solutions favored. The global anal- 
ysis of the solar and reactor neutrino results yields 
Am 2 = 71±l ° 3 x 1CT 5 eV 2 and = 32.51^ degrees. 

SNO is presently analyzing its full salt data set with 
a detailed treatment of the day /night and spectral 
information. In the future SNO will perform a global 
oscillation fit with a maximum-likelihood method. 
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6. Pull Analysis 

The pull method allows a split of the residuals 
from the observables and the systematic uncertain- 
ties 01- This alternative approach embeds the effect 
of each independent k th source of systematics through 
a shift of the difference (i?^ xp — R^ 1 ) by an amount 
£fcC^. The normalization condition for the K indepen- 
dent sources of systematic uncertainty is implemented 
through quadratic penalties in the global x 2 ? which is 
minimized with respect to all e^'s 
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In an experimental context, the pull approach is not 
blind since it uses the data to constrain the systematic 
uncertainties. Systematic shifts calculated with the 
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fit after a detailed study of the systematic uncertainty 
of each observables. See details in Ref. [l4| . 



7. Summary 
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